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ABSTRACT 



We have obtained new Tully-Fisher (TF) peculiar velocity measurements 
for 52 Abell galaxy clusters distributed throughout the sky between ~50 and 
200/i~ 1 Mpc. The measurements are based on / band photometry and optical 
rotation curves for a sample of 522 spiral galaxies, from which an accurate TF 
template relation has been constructed. Individual cluster TF relations are 
referred to the template to compute cluster peculiar motions. 

The reflex motion of the Local Group of galaxies is measured with respect 
to the reference frame defined by our cluster sample and the distant portion 
of the Giovanelli et al. (1998a) cluster set. We find the Local Group motion 
in this frame to be 565±113 km s" 1 in the direction (7,6) = (267°, 26°) ± 10° 
when peculiar velocities are weighted according to their errors. After optimizing 
the dipole calculation to sample equal volumes equally, the vector is 509 ± 195 
km s" 1 towards (255°, 33°) ± 22°. Both solutions agree, to within la or better, 
with the Local Group motion as inferred from the cosmic microwave background 
(CMB) dipole. Thus, the cluster sample as a whole moves slowly in the CMB 
reference frame, its bulk flow being at most 200 km s _1 . 
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Subject headings: galaxies: distances and redshifts — cosmology: observations; 
distance scale 



1. INTRODUCTION 

The motion of the Local Group of galaxies is imprecisely known. Using the COBE 
result for the motion of the Sun in the CMB reference frame (Kogut et al. 1993) and 
the recipe of de Vaucouleurs et al. (1976) to convert heliocentric velocities to the Local 
Group frame yields a Local Group velocity of Vi g _ cm b = 620 ± 22 km s _1 towards 
(I, b) = (271°, 29°) ±3°. The volume over which the coherent motion represented by Vi g _ cm b 
persists is unknown. It can be estimated, however, by measuring the reflex motion, with 
respect to the Local Group, of spherical shells of increasing radii (Giovanelli et al. 1998b; 
G98b). The "convergence depth" is the distance at which the shells' dipole motions begin 
to mirror the Local Group motion. Alternatively, the distance at which the shells' bulk 
motions, the dipole motions in the CMB reference frame, approach rest defines convergence. 
Using Tully-Fisher (TF) distances for a large sample of field galaxies out to 95/i" 1 Mpc, 
G98b find the convergence depth to be of order 60/i _1 Mpc. This result is strengthened by 
the work of Giovanelli et al. (1998a; G98a), who find evidence for a small bulk flow using 
TF distances on a sample of clusters between 30 and 92h~ 1 Mpc. These efforts, however, do 
not preclude non-zero bulk motions on scales as large as those claimed by Lauer & Postman 
(1994; LP). In this letter we tackle the issue of large-scale bulk motions by measuring 
the reflex motion, with respect to the Local Group, of a deep, all-sky sample of 52 Abell 
clusters that extends to ~200/i _1 Mpc. To fortify our results, we merge to that sample the 
distant half of the G98a sample. In the following section we highlight the general properties 
of the data. Section |] deals with dipole calculations of the cluster peculiar velocity sample 
and Section |] covers the effect the Zone of Avoidance has on such computations. We review 
and discuss our results in Section |. 

2. THE ABELL CLUSTER INERTIAL REFERENCE FRAME 

The TF data for this project have been collected over the last four years. Installments 
of the data have been presented previously (Dale et al. 1997, 1998), with the third and final 
installment forthcoming. A description of the entire sample's properties, the construction of 
the / band TF template, the sample biases and their corrections, and the individual cluster 
peculiar velocities will appear in a separate work. We briefly summarize the data here. 
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The sample consists of an all-sky collection of 52 rich Abell clusters between ~50 
and 200/i~ 1 Mpc (hereafter the 'SCII' sample). Figure [i] displays the cluster distribution 
in Galactic coordinates. The symbol sizes are scaled inversely proportional to the cluster 




Fig. 1. — The all-sky distribution of the SCII sample in Galactic coordinates. The circles 
are drawn inversely proportional to the cluster redshifts. The two examples in the lower left 
give the scale. Asterisks mark the apices of the motion of the Local Group with respect to 
the CMB and the LP cluster inertial frame. 

redshifts. The average of the cluster CMB velocities is 12,100 km s -1 . TF peculiar 
motions have been estimated for these clusters based on / band photometry and optical 
rotation curve work for 522 spiral galaxies. The number of galaxies with TF data and their 
dispersion determine the accuracy of a given cluster's peculiar motion estimate. For our 
sample, the average cluster peculiar velocity uncertainty is ~6% of the cluster distance. 

The sky and redshift distributions of our cluster sample play important roles in 
this study. Adequate sky coverage is necessary for any dipole measurement and the 
redshift distribution must be tailored to test claims of bulk flows on large scales. Our 
sample has been selected according to these precepts, but it is of course useful to consider 
supplementing our data with peculiar velocity measurements of other clusters. The distant 
members (cz > 4500 km s -1 ) of the G98a 'SCI' sample afford such an opportunity. Our 
data represent an extension of the SCI data to higher redshifts so we have taken care to 
process the data in a manner consistent with the data reductions for the SCI program. To 
ensure uniform sky coverage and a redshift overlap with the SCII sample, we restrict the 
following dipole computations to involve only the distant half of the SCI sample, the 12 
clusters beyond 45/i _1 Mpc. We shall refer to the combined sample of 64 SCI+SCII clusters 
as the SC sample. The average redshift of the SC sample in the CMB frame is anywhere 
from 7500 to 11,500 km s" 1 , depending on the weighting scheme. We point out that this 
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sample depth matches that of the LP data set, a sample for which the local bulk flow was 
found to be on the order of 700 km s _1 . 



The reflex motion of the Local Group with respect to the SC sample, Vi g _ sc , follows 
from a chi-squared minimization of the dipole vector with respect to the cluster sample's 
peculiar velocity field in the Local Group frame, i.e. the Vi g _ sc that minimizes 



where fj is the unit vector in the direction of the i cluster, and Vi and are 
the cluster's peculiar velocity and its error. The factor 1/si allows us to vary the 
weighting scheme; we assume l/sj=l in what immediately follows, and we will 
henceforth call this "Case a." A straightforward minimization of this equation yields 
Vig_sc = (V x , V y , V z ) = (—29, —507, 247) km s" 1 in Galactic Cartesian coordinates, or 
alternatively, 565 km s _1 towards (l,b) = (267°, 26°). A robust approach to estimating the 
uncertainty in the dipole is to utilize the characteristics of the cloud of solutions generated 
from Monte Carlo simulations. In these simulations we generate a Monte Carlo trial sample 
by adding random Gaussian deviates of the peculiar velocity errors to the original peculiar 
velocities. The dipole solution is calculated, and the process is repeated for a large number 
(10 4 ) of realizations. The dispersions in the dipole components, (93,119,80) km s _1 , imply 
an error-bias corrected dipole amplitude of 538 ±113 km s _1 (error bias systematically 
inflates the dipole amplitude: V 2 m = V 2 + V 2 + V 2 — e 2 — e 2 — ef). Figure |] displays the 
error cloud projected onto the three planes of the Galactic Cartesian coordinate system. 
We can use this error cloud to better understand the statistical likelihood domains of 
the measurement; the contours plotted enclose 68% and 95% of the points, respectively 
representing la and 2a confidence regions. The contours derive from two-dimensional 
Gaussian fits to the Monte Carlo distribution. The filled circles indicate the dipole motion of 
the Local Group with respect to the CMB, whereas the open circles with the accompanying 
error bars show the value as determined by LP. 

A Ax 2 analysis can be another approach to understanding the dipole result. The 
difference in \ 2 using our solution and one with a null reflex Local Group motion is 
Ax 2 = 51.4 (or Ax 2 = 38.5 if we account for the errors on the CMB dipole). For a three 
parameter fit this implies we can guarantee non-zero Local Group motion with 99.999% 



3. 



DIPOLES AND BULK FLOWS 



3.1. The Abell Cluster Dipole 
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Fig. 2. — Galactic Cartesian coordinate representation of the motion of the Local Group 
with respect to the Abell cluster inertial reference frame. The small dots plotted are Monte 
Carlo realizations of the peculiar velocity data set, and the ellipses drawn represent the la 
and 2a confidence regions of the dipole. The large filled circles indicate the motion of the 
Local Group with respect to the CMB and the open circles show the dipole motion with 
respect to the LP reference frame, along with the quoted uncertainty. 

confidence. More importantly, we can test for agreement with the COBE solution for Local 
Group motion, Vi g _ cm b = (6, —541, 303) km s -1 . In this case A% 2 = 1.4; our results for the 
Local Group motion agree with that of COBE to well within la — the COBE solution lies 
within our 68% confidence region. Conversely, A% 2 = 17.8 in a comparison between our 
solution and that of LP (computed including the LP error estimate); we disagree with the 
LP result to better than 99.9% (a 3.5a difference). 

Table 1 summarizes the dipole results, parenthetically including the dipole amplitudes 
corrected for error biasing. An indication of the effective sample depth is included, 
computed by averaging cluster recessional velocities weighted by l/(sje 2 ). The number of 
clusters considered in each computation, N c \, exceeds the number of SC clusters (64) when 
artificial clusters from the Zone of Avoidance are included, as described later. Finally, we 
have also computed the dipole after correcting for a homogeneous Malmquist bias, in the 
form advocated by G98b. This correction alters results little; it is calculated to estimate 
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the amplitude and potential impact of the effect. It is likely that a more complicated, 
but harder to quantify, inhomogeneous bias correction is necessary. However, the variance 
between the inhomogeneous and the homogeneous forms of the correction is likely to be a 
second order effect; either form impacts little the final result. 

Table 1. Cluster Dipole 
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"homogeneous Malmquist bias correction 

6 Artificial ZoA clusters included, with null average peculiar velocity in Local Group frame 
c Artificial ZoA clusters included, with null average peculiar velocity in LP frame 



3.2. The Dipole Under Equal Volume Weighting 

The previous formalism does not weigh equal volumes equally, but rather tends to 
weigh nearby regions more heavily. This unequal volume weighting exists because the 
completeness of our cluster sample and the accuracy of our peculiar velocities fall with 
increasing distance. An exponential fit to the completeness function of the SCI+SCII 
sample follows exp(— cs[km/s]/4700). To accomplish equal volume weighting, or more 
precisely equal spherical shell weighting, we construct the factor in Equation [1] to 
account for the characteristics of our cluster selection function. The resulting combination 
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of our cluster selection function and the peculiar velocity error distribution yields a roughly 
flat weighting with redshift: l/(s*e?) oc 1/y/Ni, hereafter referred to as "Case 6," where Ni 
is the number of galaxies used in the i th cluster. Results for the dipole computations using 
this fair weighting scheme are also listed in Table 1. We again find agreement with the 
COBE solution, though the associated errors are admittedly larger - giving similar weights 
to near and far clusters inevitably increases the dipole uncertainty since distant clusters 
have less accurate peculiar velocities. Results are also listed for the dipole computation 
when the cluster Abell 3667 is excluded, an object with a highly uncertain peculiar velocity 
estimate, as it is based on TF measurements of only four galaxies. By performing jackknife 
and statistical tests of the Case b dipole, we find this cluster to unduly alter the results (its 
presence in Case a solutions makes little difference due to its large peculiar velocity error). 
We thus conservatively exclude this cluster from further Case b results. The Monte Carlo 
realizations for the Case b dipole (Figure |3]) are in good agreement with the CMB dipole. 
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Fig. 3. — Similar to Figure |2|, but for equal volume weighting: = 1. 
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3.3. The Bulk Flow 

We use these results in concert with the estimate of the Local Group motion in the 
CMB frame to ascertain the motion of the cluster inertial frame with respect to the CMB, 

V sc _ cm b = Vlg— cmb ^lg— sc- (2) 

We find the bulk flow amplitude for the two cases 1/sj = 1 and l/(sjef) oc l/y/Nl to be 
75 ± 92 km s" 1 and 178 ± 196 km s _1 ; corrected for error biasing, the amplitudes reduce to 
a null bulk flow. The inertial frame defined by the distant SCI+SCII clusters is determined 
to be at rest with respect to the CMB fiducial rest frame. We note that the above estimates 
for the local Universe's bulk flow amplitude match well what is expected from theoretical 
considerations (e.g. Feldman and Watkins 1994; Strauss et al. 1995; Moscardini et al. 1996). 



4. THE EFFECT OF THE ZONE OF AVOIDANCE 

The region of the sky that is obscured by the plane of the Galaxy is rather large even 
in the / band; the Zone of Avoidance (ZoA) extends to about ±20° in Galactic latitude, 
or alternatively, it covers about a third of the sky. We have abstained from working at 
low Galactic latitudes because Galactic extinction corrections there are highly uncertain. 
We first turn to Monte Carlo simulations to gauge the effect of the ZoA on bulk motion 
measures. For each of 10 4 trial realizations of the dipole we generate 15-20 artificial clusters 
between —20° < b < 20° that mimic the overall characteristics of the observed cluster 
sample. Table 1 includes the average results of these simulations, with the quoted errors 
deriving from the la dispersions. The amplitudes and apices agree less well with Vi g _ cm b, 
but we emphasize that the artificial clusters added have, on average, zero peculiar velocity 
in the Local Group frame in one case or they have been assigned, on average, the LP dipole 
motion in another case. This exercise is carried out to gauge, very conservatively, the 
maximum potential effect of the unsampled ZoA region. 

Numerical simulations based on the truncated Zel'dovich approximation kindly 
provided by S. Borgani allow us to further test the impact of the ZoA and the effects of 
cosmic scatter on our dipole calculations. Such simulations allow one to gauge uncertainties 
introduced by blanking large swaths of the sky. In other words, comparisons can be made 
between dipoles that sample the entire sky and those that are limited by a ZoA. Moreover, 
simulated peculiar velocities have no "observational errors." Thus, we can test for different 
effects to a higher accuracy than would otherwise be attainable. As described in Moscardini 
et al. (1996), the simulations use the Zel'dovich approximation to mimic gravitational 
dynamics. The simulations are carried out on a 256 3 grid, one particle per grid point, within 
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a cube of 960/i 1 Mpc sides, a region large enough to contain all the waves responsible for 
motions on scales of 200h~ 1 Mpc. A variety of cosmological models are simulated. 

To enable a fair comparison with our cluster observations, we cull the number of density 
peaks to a density similar to that found in the Abell/ACO catalog: n c \ ~ 2xl0 _5 (/i _1 Mpc) -3 , 
or an average cluster-cluster spacing of ~ 40/i _1 Mpc. We also only calculate a cluster 
dipole for "Local Group-like" observers. In other words, |V[ g | = 620 ± 44 km s _1 for a 
top-hat sphere of radius 7.5/i _1 Mpc, and a matter overdensity range of —0.2 < Sp/p < 1. 
Such criteria are satisfied by 10 2 to 10 3 grid points, depending on the cosmological model. 

The dipole calculation is executed three times for each Local Group-like observer 
using a random sample of clusters whose space distribution is fashioned after the combined 
SCI+SCII cluster sample. The first computation tests only for the effect of cosmic variance 
in the peculiar velocity field and thus does not involve any sort of sky masking. The apex 
of this first dipole computation can also serve as a working reference point for the "true" 
dipole. The second calculation includes a randomly placed ZoA for each observer, whereas 
the third calculation involves a ZoA mask that is oriented, with respect to the true apex, in 
the same manner we observe the plane of the Milky Way with respect to the apex of our 
own Local Group motion in the CMB frame. The variance of the dipole solution for the 
set of Local Group-like observers gives us an idea of how well we can recover dipole flows. 
The majority of the simulations indicate that the measurement of 60-70 cluster cosmic 
motions limit bulk flow measurements to an accuracy no better than ~100-150 km s -1 ; the 
presence of a ZoA further confuses dipole measurements by ~60 km s^ 1 , resulting in an 
overall uncertainty of ~120-160 km s _1 , consistent with the errors in Table 1. 



5. CONCLUSIONS 

Our goal is to conclusively test claims of high amplitude, large scale bulk motions 
using Tully-Fisher redshift-independent distances to galaxy clusters. In collaboration with 
the accurate SCI peculiar velocities, we use our sample of 52 cluster peculiar velocities to 
accomplish our goal, and in the process we establish a reliable upper limit to the local 
convergence depth. We briefly recap here our findings. 

• The reflex motion of the Local Group is calculated with respect to the Abell cluster 
inertial reference frame represented by the distant SCI+SCII sample. Upon weighting 
peculiar velocities by their error estimates, we determine the Local Group motion to be 
565±113 km s" 1 in the direction (Z,6) = (267°, 26°) ± 10°. Employing a different weighting 
scheme, one where equal volumes are treated equally, yields 509 ± 195 km s" 1 towards 
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(255°, 33°) ± 22°. Both results agree with the CMB dipole. The bulk motion of the SC 
sample in the CMB reference frame is at most 200 km s _1 . 

• We employ a variety of numerical simulations to test the statistical validity of our 
results. We find that bulk flow measures not derived from a truly all-sky sample, but rather 
those from samples suffering from a dearth of objects at low galactic latitudes such as for 
our survey, can be reliable. This fact holds for a variety of cosmological scenarios. 

• We have applied a suite of statistical tests to our dipole calculations. Bootstrap and 
subset resamplings in addition to jackknife tests indicate our results are robust. 

• We find evidence supporting the small convergence depth of order 60h~ l Mpc found 
by Giovanelli and coworkers. They use two independent samples of TF distances, one 
comprised of 24 clusters out to 92/i~ 1 Mpc (G98a) and the other made up of a large sample 
of spiral field galaxies covering a similar distance regime (the 'SFF sample; G98b). With 
respect to spherical shells of increasing radii, they find the reflex motion of the Local Group 
increasingly mirrors the CMB dipole. Our result shows that the bulk flow of the volume 
subtended by the SC sample, of effective depth ~100/i _1 Mpc, is less than 200 km s" 1 , thus 
confirming convergence is maintained beyond the limits of the SCI and SFI samples. 
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